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O ', Abstract: This paper focuses on the analysis of the gravitational instability in 

presence of bulk viscosity both in Newtonian regime and in the fully-relativistic 
approach. The standard Jeans Mechanism and the Quasi- Isotropic Solution are 
treated expressing the bulk-viscosity coefficient £ as a power-law of the fluid 
energy density p, i.e., C — CoP s - In the Newtonian regime, the perturbation 
evolution is founded to be damped by viscosity and the top-down mechanism 
of structure fragmentation is suppressed. The value of the Jeans Mass remains 
unchanged also in presence of viscosity. In the relativistic approach, we get a 
\q \ power-law solution existing only in correspondence to a restricted domain of £ . 
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1 Characterization of viscosity 



The aim of this work is to study the dynamics of the gravitational instability i.e., the 
evolution of scalar density perturbations, in presence of viscous effects. Viscosity can be 
divided in two different kinds: bulk viscosity (Bv ) ( and shear viscosity In this paper, we 
will concentrate on the analysis of isotropic (or almost quasi-isotropic) cosmological models. 
In this respect, we can safely neglect shear viscosity since there is no displacement of matter 
layers wrt each other (in the zeroth-order motion) and this kind of viscosity represents the 
energy dissipation due to this feature. Indeed, in presence of small inhomogeneities, such 
effects should be taken into account but we focus on the study of scalar density perturbations 
and volume changes of a given mass scale are essentially involved. We expect that the non- 
equilibrium dynamics of matter compression and rarefaction is more relevant than frictions, 
thus, we concentrate on Bv only since it outcomes from the difficulty for a thermodynamical 
system to follow the equilibrium configuration. 

In this paper, we assume Bv as function of the state parameters of the fluid following the 
line of the fundamental analysis due to the Landau School pQ. In particular, we implement 
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the so-called hydrodynamical description of the fluid, i.e., an arbitrary state is consistently 
characterized by the particle-flow vector and the energy momentum tensor alone and viscosity 
is fixed by the macroscopic fluid parameters. In the homogeneous models £ depends only on 
time and we assume it as a power-law of the density of the fluid p, i.e., 

C = CoP s , (i) 

with (o, s = conts. 



2 Jeans Instability in Presence of Bulk Viscosity 

The study of the density-perturbation dynamics during the matter-dominated era can be 
consistently described using the Newtonian picture, as soon as sub-horizon-sized scales are 
treated. The fundamental result is the Jeans Mass [5] (Mj), which is the threshold value 
for the fluctuation masses to condense, generating a real structure. If masses greater than 
Mj are addressed, the density contrast (6) diverges (in time) giving rise to the gravitational 
collapse. 

In the following, we brief discuss how Bv affects such a scheme. The starting point 
is the Eulerian set of motion equations on which a perturbative theory is developed (the 
background is assumed static and uniform). If we introduce Bv effects in the first-order 
analysis, the density-contrast growth results to be dumped by viscosity, suppressing the 
structure formation without changing the threshold value of Mj. In particular, the grater £ 
is, the slower 5 diverges in time and, furthermore, viscosity generates a decreasing exponential 
regime and a damped oscillatory one in place of the standard pure oscillatory behavior [3] . 

In this respect, we now analyze the top-down mechanism of sub-structure formation, i.e., 
the comparison of one collapsing agglomerate with M 3> Mj and an internal non-collapsing 
sub-structure with M < Mj. The sub-structure mass must be compared with a decreasing 
Mj since the latter is inversely proportional to the collapsing agglomerate background mass. 
As soon as such value reaches the sub-structure mass, this begins to condense. In the stan- 
dard Jeans Model, such mechanism is always allowed since the amplitude for perturbations 
characterized by M < Mj remains constant in time but, the presence of decreasing S in 
the viscous model, requires a discussion on the effective damping and on the efficacy of the 
top-down scheme. As shown in Figure 1, in correspondence of a very small viscosity coeffi- 
cient (case A), we can show how the sub-structure survives in the oscillatory regime during 
the background collapse and the fragmentation occurs. On the other hand, if Bv is strong 
enough (case B), the damping becomes very strong and the sub-structure vanishes during 
the agglomerate evolution. Thus, the top-down mechanism results to be deeply unfavored. 
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Figure 1 : A - Co 

B - Co ■ 



- 10 g/cms, galaxy 10 Mq: Sq (dashed), sub-structure 10 Mq: 8s (normal). 
14 g/cm s, galaxy 10 12 Mq : 8q (dashed) , sub-structure 1 Mq : 5s (normal). 



Moreover, if the static and uniform background solution is corrected for the expansion of 
the Universe [3J, a Jeans- like relation and a considerable damping of the density contrast 
growth can be found. 
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3 Bulk-Viscosity Effects on the Quasi-lsotropic Solution 



In 1963 [3], E.M. Lifshitz and I.M. Khalatnikov firstly proposed the Quasi-lsotropic Solution 
(QIS) based on the idea that, as a function of time, the 5-metric (where a, /3=1, 2, 3) 
is expandable in powers of t, i.e., a Taylor expansion is addressed: 

7a/3 = t x a aP (x) + i» b a0 (x) , 7 a/3 = t~ x a a ? - , (2) 

where x > (space contraction) and y > x (consistence of the perturbation scheme) . 

We now focus on the relevance of dealing with Bv properties of the cosmological fluid 
approaching the Big-Bang singularity 5 . As far as we characterize the Bv coefficient like 
C = ( p s , it is easy to realize that the choice s — 1/2 prevents dominating viscous effects. 
On the other hand, simple considerations indicate that the case s < l /a leads to negligible 
contributions of Bv in the asymptotic regime towards the Big-Bang. As a consequence, in 
studying the singularity physics, the most appropriate form of the power-law is C = Co \[P an d 
our aim is to determine the conditions on the parameter Co (i.e., on the viscosity intensity) 
which allows for the existence of a QIS for the radiation dominated Universe. 

To this purpose, we separate zeroth- and first-order terms into the 3-metric tensor and the 
whole analysis of Einstein equations follows this scheme of approximation. In the search for 
a self-consistent solution, we make use of the hydrodynamical equations in view of fixing the 
form of the energy density. Of course, in our solution the power-law for the leading 3-metric 
term is sensitive to Co- As a result, by guaranteeing the consistence of the model, we find 
that the QIS exists if and only if Co remains smaller than a certain critical value, i.e., 

Co < Co* = 2 /3V3 , (3) 

(here 8irG = c = 1). In fact, for values of the viscous parameter greater than Co 5 the 
perturbative expansion towards the singularity can not be addressed, since fluctuations would 
grow more rapidly than zeroth-order terms. 

It is worth noting that the study of the perturbation dynamics in a pure isotropic picture 
yields a very similar asymptotic behavior when viscous effects are taken into account [B] . The 
Friedmann scheme is preserved only if we deal with limited values of the viscosity parameter, 
obtaining the condition ^ lso ^ < Q /3. This smaller constraint is physically motivated if we 
consider, as it is, the Friedmann model a particular case of the QIS. Finally, analyzing S, we 
confirm and generalize the result obtained in [B] about the damping of density perturbations 
by the viscous correction. 
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